Magnetic Catalysis in the Presence of Scalar Fields 
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The influence of scalar fleld interactions in the magnetic catalysis phenomenom is studied in a 
gauge theory with scalar and fermion fields. It is shown that the external magnetic field catalyzes the 
appearance of a fermion-antifermion condensate along with a non-zero scalar vev, hence generating 
a fermion dynamical mass and breaking the discrete chiral symmetry of the theory. The model does 
not require the introduction of a scalar mass term with a wrong sign in the original Lagrangian to 
generate the scalar vev, nor does it exhibit dimensional transmutation "a la Coleman- Weinberg." 
. . . The scalar vev noticeably enhances the dynamically generated fermion mass. Possible cosmological 

^— «' , applications of these results are suggested. 

PACS: ll.SOQc, 11.30 Rd, 12.15-y. Keywords: non-perturbative effects, external magnetic fields, 
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Symmetry behavior in quantum field theories under the influence of external fields has long been a topic of intensive 
study in theoretical physics . In the present paper we are particularly interested in non-perturbative effects produced 
, by external magnetic fields in gauge theories with scalar-scalar and scalar- fermion interactions. 
^< ■ When a non-perturbative analytic approach is used to study field theories in external magnetic fields, new non- 
CNJ ' trivial effects are found. An important example of these non-perturbative effects is the formation of a chiral symmetry 
i breaking fermion condensate < ipip > and of a dynamically generated fermion mass in the presence of an external 
magnetic field, known in the literature as magnetic catalysis Q, As was originally shown by Gusynin, Miransky 
and Shovkovy |^ , the magnetic field catalyzes the symmetry breaking by reducing to the weakest attractive coupling 
' the strength of the interaction needed to break the symmetry. This phenomenon has proven to be rather universal 
, and model independent The magnetic catalysis is not only interesting from a purely fundamental point of 

^—1 ■ view, but it has potential application in condensed matter and cosmology ^. 

In the present paper we aim to investigate the infiuence of scalar field interactions in the magnetic catalysis 
I phenomenon. With this goal, we will consider a simple model field theory containing self interacting scalar fields, 
I as well as Yukawa fermion-scalar interactions in the presence of an external constant magnetic field. Using a self- 
. consistent approach, we will determine how the existence of these interactions in the presence of the magnetic field can 
' affect the stability of the vacuum. We will see below that due to the magnetic field, a non-zero vacuum expectation 
Oh, value (vev) of the scalar field and a fermion dynamical mass arise as the solution of the minimum equations of the 
f-^ • system, breaking in this way the discrete chiral symmetry of the original Lagrangian. However, in contrast to more 
' conventional mechanisms to generate scalar vev's, the present model does not require the introduction of a scalar 
■ mass term with a wrong sign in the original Lagrangian, nor does it have dimensional transmutation " a la Coleman- 
^ I Weinberg." Instead, the scalar vev is catalyzed, along with a fermion-antifermion condensate, by the external magnetic 
• ^ 1 field. The only assumption we make is that all couplings are weak enough to justify the Hartree-Fock approach used 
' below. 

Both the fermionic condensate and the scalar vev contribute to the dynamically generated fermion mass. We find 
' that the scalar vev and the fermion dynamical mass grow with the square root of the magnetic field strength. The 
main effect of the scalar interactions is a significant enhancement of the magnetic catalysis as reflected in the increase 
of the fermion dynamical mass compared to the case without scalars. The mass enhancement resembles a similar one 
produced by lowering the number of spatial dimensions in the Nambu-Jona-Lasinio (NJL) model from 3-1-1 to 2-1-1 
Like other theories with magnetic catalysis, no critical magnetic field is required in the present model for the 
zero-temperature chiral symmetry breaking to take place. 

The studied model is not intended as a realistic theory, but rather as an example of a large class of theories with 
scalar fields, on which a dynamical symmetry breaking can be catalyzed by an external magnetic field. In this sense, 
it could be useful for condensed matter, as well as for cosmological applications. 

Our starting point will be a theory of gauge, fermion and real scalar fields described by the following Lagrangian 
density 

L = - jF-'F^, + iin'Sf* + 9fn''*K - \s^vS''v - j*)' - y*)' - A.vJ** (1) 
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It has a U(l) gauge symmetry, 



a fermion number global symmetry 



and a discrete chiral symmetry 



^ ^ e*"(^V, (2) 



- e^V, (3) 



V'^TbV', ^'^-'075- ip-^-ip (4) 



Note that a fermion mass term mijj^p is forbidden, since it is invariant under (g) and (||), but not under the discrete 
chiral symmetry (^. 

To study the vacuum solutions that could arise in the theory under the influence of an external constant 
magnetic field B (without loss of generality we will assume throughout the paper that sgn{gB) > 0), we need to solve 
the extremum equations of the effective action F for composite operators 

%^^0, (5, 

oG 

'li^^O (6) 

where G{x,x) — (t{x) = ^0 | ■4'{x)4'{x) | O) is a composite fermion-antifermion field, and ipc represents the vev of the 
scalar field. If the minimum solutions of (|^) and (^) are different from zero the discrete chiral symmetry (Q) will be 
broken and the fcrmions will acquire mass. We shall see that due to the magnetic field this is indeed the case. 

We choose the composite field G{x,x), ignoring the second possible one, n{x) — (O | ip{x)i"f^'4'{x) \ O) , since the 
effective action can be a function only of the chirally invariant combination = cr^ + tt^. 

The loop expansion of the effective action F for composite operators |jll[, can be expressed as 

F (G, ipc) = S [ipc) - iTr In G^^ + i]^Tr In + i^Tr In A^^ - iTr [G^^ {ipc) G] + F2 (G, ^c) + C (7) 

Here G is a constant and S (ipc) is the classical action evaluated in the scalar vev ipc- The bar on the fermion propagator 
G {x, y) means that it is taken full, while the non-bar notation indicates free propagators, as it is the case for the gauge 

propagator D^y{x -y)^ J '"'^l^I ' (sm"- - (1 - O^zij) > and the scalar one A{x -y)^ J -^-$^^M^e^ ^it^^ 
A'P — ^(p'^ + is^. A dependence on full boson propagators is not included, since we do not expect that any of the boson 
fields will acquire expectation values of their composite operators. In general F2 (G, (pc) represents the sum of two 
and higher loop two-particle irreducible vacuum diagrams with respect to fermion lines. For a weakly coupled theory 
one can use the Hartree-Fock approximation, which corresponds to retaining only that contribution to F2 which is 
lowest-order in coupling constants (i.e. two- loop graphs only). In this case F2 is 



2 /> 

r2(G,^c) = Yj d^xd^ytr\G{x,y)-f^G{y,x)YD^,{x,y) 



I d^xd^ytr (x, x)) D^^x - y)tr {^''G (y, y)) 
I d^xdStr[G{x,y)Giy,x)A{x,y)] 

\2 

2 



d^xd^ytr (G (x, x)) A(a; - y)tr (G (y, y)) (8) 



The extremum equations (|^) and (^) correspond, respectively, to the Schwinger-Dyson (SD) equation for the fermion 
self-energy operator E (gap equation), and to the usual minimum equation for the expectation value of the scalar 
field, which in the presence of the magnetic field has to be determined in a self-consistent way, that is, simultaneously 
with the gap equation. 
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Although we have introduced a bare scalar mass //in (|^), as we are interested in the possibility of a dynamically 
generated scalar vev, we shall take the limit // — > at the end of our calculations. 

The second to fifth terms in the effective action (Q) correspond to one-loop contribution. Their evaluation is quite 
straightforward (the scalar self- interaction can be renormalized in the usual way), with the exception of the fcrmion 
contributions, which contain the background magnetic field. Then, let us calculate explicitly the one-loop fermion 
contribution coming from the term 

r« = -^rr[G-i(^,)G] (9) 

in (|7|). The trace operator is along continuous and discrete variables. 

The full fermion propagator G{x,y) in the presence of a constant magnetic field B appearing in Eq. (|^) can be 
written as Hi, [| 

Similarly, the free fermion inverse propagator in the presence of _B, denoted by (ipc) , is given by 

G-i(^,y,V.c)=E / ^4?^^P'(^)(7-P +™o)^,,(y), (11) 

with p ~ (pojO, ~\/2gBk,p3), k the Landau level index, and mo — Xy(pc the fermion mass appearing after the shift 
tf ^ If + (pc in the scalar field. The magnetic field B has been taken along the third axis. 

In the above equations we have introduced the Ritus' Ep functions ||l^. These orthonormal and complete set of 
function-matrices provides an alternative method to the Schwinger's approach to problems of QFT on electromagnetic 
backgrounds^ Details on the Ep functions properties and the notation used below can be found in ref. 

Using Eqs. ( p^ and (pT|), and the orthonormality of the Ep functions, the integrals in the coordinate variables, 
as well as the sum and integrals in one of the momentum variables can be done straightforwardly to find that the 
one-loop fermion contribution (O) can be expressed as 



(2^) 

where the notation 5'^^^(fc) = 6{ko)S{k2)6{k:i) is understood. 

At this point we need to consider the structure of the mass operator E introduced in ref. 

J2ip) = (Ph ■ + Z±{p)j ■ p^ + m{p) (13) 

where p^^ — {pqO, 0,^3), p^ — (0, 0, —y/2gBk, 0). The coefficients Z^^ (p), Z±{p) and m(p) are functions of p^. Note the 
usual separation in the presence of a magnetic field between parallel and transverse variables. Substituting with the 
structure ( |l3|) in Eq. (p^, taking the trace, and performing the Wick rotation to Euclidean coordinates, we obtain, 
after integrating in p2 

rW 8 RA(4)rm\- f, , il + Z,)pl + {l + Z^)pl+m{p)mo 

We are interested in the contribution of T^p to the minimum equations (||) and (^). In the case of the gap equation, 

such a contribution can be found differentiating directly Eq. (^ with respect to G. For the minimum of the scalar 
field, we need the derivative of (|l^) with respect to ipc- It takes the form 

= 87TXygB5(^\0) J2 / dp,dp, _ pP>_ (15) 



^For an application of Ritus' method to the QED Schwinger-Dyson equation in a magnetic field see ref. |Q. 
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where the solutionis Z^^ = of the SD equation (||) was explicitly used. 

For large magnetic field, the main contribution to the sum in k comes from the lower Landau level, i.e. = 
2gBk = 0. Within this approximation Eq. ( p5| ) becomes 

= 8nXygBS^^\0) / dp,dps _^ ' l'_ (16) 
0(pc J p|+m^(p[|) 

In general, the dynamical mass depends on the momentum. Nevertheless, it is reasonable to expect that, 

similarly to QED ||^, m{pii) behaves as a constant in the infrared region, and diminishes with increasing 



< \/gB. Therefore, we can 



Then the integral (|l^) is dominated by the contributions from the infrared region 

approximate the function r7i(p||) by a constant solution to(]3||) ~ m(o) — m <^ \/gB and cut the momentum integration 
at ^/gB to obtain 

^ = V(^)^gBm\n ( 4) - -V^''' Xy < > (17) 

where V'^^^ represents an infinite four dimensional volume, and < -0?/; >= irr|G(x,x)} — — In ^f^^ denotes 

the fermion-antifermion condensate ^ induced by the external magnetic field. 

The two-loop contributions are a little more involved. As we have not enough space in a letter to give all the 
detailed calculations, we will explicitly show, for the sake of understanding, the evaluation of one term. The others 
can be found in a similar way. The complete calculation will be published elsewhere. 

First, notice that the second and fourth term in Eq. (^ generate tadpole diagrams in the SD equation (^). It is 
easy to realize that the tadpole diagram with the gauge-fermion vertex vanishes. However, the tadpole associated to 
the scalar-fermion vertex is not zero when (^c 7^ and, as shown below, it has a significant contribution to the gap 
equation. Let us evaluate this tadpole contribution, which we denote by ■ 

H^i^, V) = = -i>?J\x - y) / dSA(x - z)tr G[z, z) (18) 



5G y 

We can transform Eq. (p^) to momentum space with the help of the Ep (x) functions to obtain 



j d^xd^yEpix) jf{x,y)Ep, (y) = (27r)4j(4) (p _ p') ^^( 
d^q e*«(^-^) 



'iXf. / d'^xd'^zEp(x) I —5 : — 



X 

k' 



and then integrate in x and z to find 

(2^)4?W(p-p')E''® - -*^'/ / dVS^'Hq)S^'\p' + q-p) 



_2 e*9i(P2+P2-2p2")/2gB gi(n-n')v 



"^"^^^Tm^TT^ 4^ Vn(fc,a)!.-(fc-,.)! ^""'^^"^^^"^ 



^The demostration that Z^^ = is a solution of the gap equation in the present theory can be done along the same line of 
reasoning followed in the Appendix of the first paper of ref. M. 
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A (a") 



7-r + E(r) 



(20) 



where = (q^ + (/f) /.gS, (ys = arctan (^^^ , and cr) = fc + | — i. 

This result can be further simpHfied by approximating the J— functions by their smaU-g^ hmit 

min(n,n") 

TO!(rt — m)!(n — mj! pi — n |! 



which is justified by the presence of the exponential factor e^'?-'- in the integrand of Eq.(|20|). Now the sums in the 
a's and the trace can be easily done, along with all but two momentum integrals, to find the following expression in 
Euclidean space, 

^ = 2^ApE / ^?'4rfP3- -2 (22) 

k" (1 + Z ) p "2 + (l + Zi)'|5,"2 + m2(p") 

V II / II 



Just like in the one-loop case, it can be seen that the main contribution to Eq. (|22| ) comes from the k" = term 
of the sum. Thus, the tadpole contribution to the gap equation reduces to 

where the solution Z — has been used. 

II 

One can recognize here the same integral that led us to the fermion condensate in Eq. iRw. Therefore, we can 
write 

_ j_^hLgBm\ii ( —\ - -2^ 

TT^ A(^2 J Xip^ 



E -^T^5S"^lnM4 =-2t^<V'V'>, (24) 



Note that the tadpole term is proportional to the magnetic field and inversely proportional to the scalar mass. We 
shall see below that such a functional dependence will be responsible for a visible increment of the mass solution in 
this model as compared to other theories. 

Taking into account all the leading contributions to Eqs. (H) and @ at large magnetic field, one arrives at the 
following minimum equations for the fermion mass and the scalar vev respectively. 



Eq. ( p^ ) can be simplified by noting that if we assume that ^ \/gB, one can neglect the term proportional 
to coming from the one-loop scalar self-interaction, compared to the term coming from the fermion condensate 
contribution ~ < >. Then, the scalar minimum satisfies 



3 Aj; 3gB ( gB 



fc^^^rn\n[^) (27) 

Likewise, the contributions of the bubble diagrams to the gap equation, (second term in (p5|)), are small compared 
to the tadpole contribution, (third term in (Eq)). Hence, equation (pa) becomes 



mo + ^-^gB7n\n[^] (28) 



Substituting with Eq. (Ej in (M) , we find 
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where the coefficient k satisfies 



-^VgB (29) 



KlnK~1.4^ (30) 



The corresponding solution for the scalar vev is 



^c^Z^V^ (31) 



Notice that the solutions ( p9| ) and (^ are indeed non-perturbative in the couplings constants. At each fixed scalar 
self-coupling A, the values of m and ipc increase with A^, since the parameter k falls down much more rapid than j- 



as Xy increases. 



It is a well known fact that in the absence of a magnetic field, the one-loop effective action (effective potential) of 
the present model would have a minimum at some non-trivial value of the scalar vev, but this minimum would lie 
far outside the expected range of validity of the one-loop approximation, even for arbitrarily small coupling constant, 
so it would have to be rejected as an artifact of the used approximation |14|. As we have just discussed, when a 



dr 



1) 



magnetic field is present, the term , Eq^(0), being proportional to the fermionic condensate, dominates the 
non-perturbative radiative corrections in Eq. (|26|). Accordingly, a non zero scalar field minimum ( [3l| ) , which, as can 
be explicitly checked out, is in agreement with the used approximation tpc ^ VgB, arises. In other words, thanks to 
the magnetic field, a consistent scalar field minimum solution is generated by non-perturbative radiative corrections. 
In this sense, a sort of non-perturbative Coleman- Weinberg mechanism takes place, with the difference that in the 
present case no dimensional transmutation occurs. Since the theory already contains a dimensional parameter: the 
magnetic field i?, there is no need to include scalar-gauge interactions in order to trade a dimensionless coupling for 
the dimensional parameter ipc [0- No constraint between the couplings has to be assumed, except that they are all 
sufficiently weak to justify the used Hartree-Fock approximation for the effective action. 

It is clear from the obtained results that the scalar field interactions significantly enhance the magnetic catalysis. 
This is better understood comparing the dynamical mass found in this work with masses generated by magnetic 
catalysis in models without scalars. For instance, in QED one has m ~ y/gB exp{— ^ ^J^) , which is much smaller 
mass than ( p9| ) (for Ay = 0.7 and A = 0.4, which correspond to the top quark Yukawa and Higgs mass Mh = llSGew, 
the dynamical mass in the theory with scalar fields is nine orders of magnitude larger). The enhancement of the 
dynamical mass due to the scalar field interactions is comparable to the effect produced by lowering the number 
of spatial dimensions in a theory that already exhibits magnetic catalysis in 3-f 1 dimensions. For example, in 

the NJL model the magnetically catalyzed dynamical mass is visibly increased from = -1^^ exp ^— ^eB^ ) ' ^^^h 

K = — in the (3-f l)-dimensional case to = \eB\^ (2-f l)-dimensional one 

A noteworthy feature of our results is that there is no trivial solution (stable or unstable) for the scalar vev in the 
present theory. Besides, no critical value of the magnetic field is required to produce the fermionic condensate and 
the scalar vev. 

When studying this model at finite temperature, it is logical to expect that the increase in the dynamical mass due 
to the scalar field interactions will lead to a corresponding increase in the temperature at which the discrete symmetry 
is restored. As known, the critical temperature separating chiral broken-unbroken phases is typically of the order of 
the zero-temperature dynamical mass. 

We can advance, although it is more or less obvious, that if the current model is extended to include a complex 
scalar instead of a real one, none of the main results will be qualitatively changed. Such a model, however, would 
resemble a simplified version of the electromagnetic sector of the SU{2) x U{1) electroweak model (the scalar would 
not be coupled to the electromagnetic field, just as the Higgs field of the electroweak theory). Since no tree- level scalar 
mass term with the wrong sign would be introduced in such a theory, the external magnetic field would catalyze gauge 
symmetry breaking by producing a non-trivial scalar vev through the non-perturbative mechanism reported in this 
paper. 

It is still too premature to envision all possible practical implications of the present work. It seems to us however 
that a careful investigation in the context of the electroweak theory of the non-perturbative effects here discussed is 
required. In the high temperature region of the electroweak theory, where the temperature-dependent fermion masses 
are nearly zero, these magnetic field-driven non-perturbative effects may be important at reasonable large magnetic 
fields. An interesting point to consider would be whether the new magnetic field dependence of the scalar vev would 
change the recent conclusions on the role of magnetic fields in electroweak baryogenesis . 
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